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Abstract — This  paper  presents  a finite  element  algorithm 
for  nonlinear  dynamic  analysis  of  cracked  beams  on  an 
elastic  foundation  subjected  to  moving  mass.  Quantity 
surveying  with  parameters  of  varied  cracks,  foundation 
and  loads  shows  their  influence  levels  on  the  nonlinear 
dynamic  response  of  the  beams.  The  findings  of  the  paper 
are  the  basis  for  the  analysis,  evaluation,  and  diagnosis 
of  damages  of  beam  structures  on  the  elastic  foundation 
subjected  to  moving  loads,  in  which  the  common  defects 
of  the  beams  such  as  cracks  are  considered  in  order  to 
improve  the  system’s  operational  efficiency  in  a wide 
range  of  engineering  applications. 
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I.  INTRODUCTION 

Beams  on  the  foundations  are  usually  modeled  to 
calculate  the  structures  of  railway  works  and  civil 
engineering.  During  the  use,  there  are  many  different 
causes  that  can  cause  weakened  defects  for  beams,  one  of 
which  is  cracks.  The  appearance  of  cracks  will  reduce 
bearing  capacity  of  the  beams,  which  leads  to  the  risk  of 
damage  to  the  building.  Salih  N Akour  [1]  analyzed  the 
nonlinear  dynamics  of  beams  on  the  elastic  foundation 
subjected  to  evenly  distributed  moving  force  by  analytical 
methods.  Also  using  analytical  methods,  Oni  and 
Awodola  [2],  Tiwari  and  Kuppa  [3]  analyzed  the 
dynamics  of  Bernoulli  - Euler  beams  on  the  elastic 
foundation  subjected  to  moving  masses.  Haitao  Yu  and 
Yong  Yuan  [4]  have  focused  on  the  analytical  solution  of 
an  infinite  Euler-Bernoulli  beam  on  a viscoelastic 
foundation  subjected  to  arbitrary  dynamic  loads.  §eref 
Dogu§can  AKBA§  [5]  investigated  the  free  Vibration  and 
Bending  of  Functionally  Graded  Beams  on  Winkler's 
elastic  foundation  using  Navier  method.  Nguyen  Dinh 
Kien,  Tran  Thi  Thom  [6]  studied  the  influences  of 
dynamic  moving  forces  on  the  Functionally  Graded 
Porous-Nonuniform  beams.  D.  Froiol,  R.  Moiolil,  E. 
Rizzi  [7]  and  D.  T.  Pham,  P.  H.  Hoang  and  T.  P.  Nguyen 
[8]  used  the  nonlinear  elastic  foundation  and  New  Non- 
Uniform  Dynamic  Foundation  applied  to  analyzed 
response  of  beam  subjected  to  moving  load  and  the 
results  show  that  the  influence  of  velocity  has  effects 
significantly  on  dynamic  response  of  structures.  N.  T. 


Khiem,  P.  T.  Hang  [9]  used  a spectral  method  applied  to 
analyzed  response  of  a multiple  Cracked  Beam  subjected 
to  moving  load. 

Using  analytical  and  finite  element  methods,  Murat.  R 
and  Yasar.  P [10],  Mihir  Kumar  Sutar  [11],  Animesh  C. 
and  Tanuja  S.  V [12],  Shakti  P Jena,  Dayal  R Parhi,  P C 
Jena  [13],  A.C.Neves,  F.M.F.  Simoes,  A.Pinto  da  Costa 
[14],  Hui  Ma  et  al.  [15]  analyzed  the  dynamics  of  cracked 
beams  subjected  to  moving  mass. 

Arash  Khassetarash,  Reza  Hassannejad  [16]  investigated 
the  Energy  dissipation  caused  by  fatigue  crack  in  beam- 
like cracked  structures.  Erasmo  Viola,  Alessandro 
Marzani,  Nicholas  Fantuzzi  [17]  used  finite  element 
method  applied  to  studied  effect  of  cracks  on  flutter  and 
divergence  instabilities  of  cracked  beams  under 
subtangential  forces. 

M Attar  et  al.  [18]  analyzed  the  dynamics  of  cracked 
beams  on  the  elastic  foundation  subjected  to  moving 
harmonic  loads  by  analytic  method,  on  the  basis  of  using 
Timoshenko  beam  model. 

So  far,  there  are  various  researches  of  beams  on  elastic 
foundation  under  transfer  (mass,  force,  oscillation 
system).  However,  for  cracked  beam  on  the  elastic 
foundation  under  moving  loads(or  masses),  most  methods 
reply  on  analytical  approaches  which  are  really  applied  to 
simple  loading  conditions.  In  this  paper,  we  develop  a 
numerical  approach  based  on  finite  element  method  for 
analyzing  the  dynamics  of  beams  on  elastic  foundation 
under  moving  masses.  We  investigate  the  influence  of  the 
elastic  foundation,  load  speeds  and  location  cracks  in  the 
dynamic  response  of  the  beams.  Note  that  finding 
analytical  solutions  of  such  beam  problems  under 
arbitrary  loading  conditions  are  really  challenging  and  no 
research  is  sufficiently  carried  out  yet.  Such  a problem 
will  be  addressed  in  this  paper. 

II.  FINITE  ELEMENT  FORMULATION  AND 
THE  GOVERNING  EQUATIONS 

A damaged  beam  has  an  open  crack  located  at  the  mid- 
section of  the  beam  at  position  x = xO.  The  beam  on  an 
elastic  foundation  described  by  an  elastic  spring  system  to 
one  direction  perpendicular  to  the  axis  of  the  beam,  which 
has  the  stiffness  kt  subjected  to  traversing  mass  4m’  at 
speed  V as  in  Fig.  1.  The  dimensions  of  the  damaged 
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beam  are  as  follows,  width  = b,  thickness  = h,  length  = L, 
crack  depth  = d. 


Fig.  1.  Cracked  beam  on  the  foundation  subjected  to 
moving  mass 
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2.1.  Equation  of  motion  of  beam  element  with  out 
crack 

The  equation  of  motion  is  derived  by  the  principle  of 
virtual  work  [19],  [21]: 

SWv+SWin+SWE=  0,  (6) 

where^Jfyis  the  virtual  work  of  internal  forces,  SWin  is 
the  virtual  work  of  inertia  forces  and  SWE  is  the  vertual 


For  finite  element  model  formulation  the  following 
assumptions  are  made:  Elastic  isotropic  materials,  cracks 
do  not  develop,  and  mass  ‘m’  is  always  in  contact  the 
surface  of  the  beam. 

The  Timoshenko  theory  describes  the  displacement  field 
components  (u,v,w)  of  an  arbitrary  point  (x,y,z)  on  the 
beam  cross-section  can  be  expressed  as 
u — Uq  (x,t^-y(pz  (x,^), 

<v  = v0(x,t),  (1) 

w = 0, 

where  uO,  vO  are  respectively  the  x and  y components  of 
the  total  displacement  vector  of  the  point  (x,0,0)  on  the 
beam  neutral  axis  at  time  t,  and  cpz  is  the  cross-section 
rotation  about  the  z-axis.  The  subscript  “0”  represents 
axis  x (y  = 0,  z = 0;  x contains  the  cross  section  centroids 
of  the  undeformed  beam,  that  will  be  often  designed  as 
middle  line  or  reference  line,  in  bending  it  coincides  with 
neutral  line).  The  x-coordinate  is  defined  along  the  beam 
length,  y-coordinate  is  along  the  height  and  the  z- 
coordinate  is  along  the  width.  The  strain-displacement 
relations  are  as 


£ = du  | 1 T dv V = dug  d<pz  y(  dvp  A 2 
* dx  2 y dx J dx  dx  2 y dx 

du  dv  dvn 

.^  = *+&=ar*’ 


(2) 


(3) 


where  |^L|  is  the  linear  part  of  the  strain  and  |^L|  is  the 


nonlinear  part  given  by: 
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The  stresses  are  related  to  the  strain  by  Hooke’s  law: 
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work  of  external  forces  due  to  a virtual  displacement. 
They  are  defined  as: 


SWV  = -J  8{s}T  {ct}  = -J  s\sf  [D]{f  ■ LjdV  - 
V V 

-j  s[eLj  [D]{eNL}dV  - J s{eNf  [£>]|^L|  dV 

V V 

-f{sNLj[D]{sNL)dV  = 


= -{Sqf  [K1]jgej  - {Sq}T  f2  ({?e})]{/}  - 


(7) 


In  this  equation  [Kl]  is  a matrix  of  constant  terms, 
[K2({q})]  and  [K3({q})]  are  matrices  that  depend  linearly 
on  the  generalized  displacements  and  [K4({q})]  is  a 
matrix  that  depends  quadratically  on  the  generalized 
displacements,  {qe}  is  the  displacements  vector.  The 
linear  stiffness  matrix  [Kl],  nonlinear  stiffness  matrices 
[K2],  [K3]  and  [K4]  have  the  following  form: 
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N= 
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(8) 


(9) 


(10) 


(5) 

where  E is  the  Young’s  modulus  of  the  material,  G is  the 
shear  modulus  and  [D]  is  the  material  matrix. 
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where  [Nu],  [Nv]  and  [N92]  are  the  row  vectors  of 
longitudinal,  transverse  along  y and  rotational  about  z 

shape  functions,  respectively,  < qe  > is  the  acceleration 


vector,  p is  the  density  of  the  beam,  [Me]  is  the  mass 
matrix,  the  vector  of  virtual  displacement  components 
will  be  represented  by  { 5d}  and  can  be  written  as 

{5d}  = {5u  5v  0}T  and  {8do}  is  the  vertual  displacements 

on  the  middle  line,  {Fo}  is  the  external  forces  on  the 
middle  line,  {Fe}  is  the  generalized  external  forces. 

The  mass  matrix  [M]  and  vector  of  generalized  external 
forces  {F}  have  following  form: 

"[M]n  [0]  [0]  ' 

[o]  [m]22  [0] 

[o]  [o]  [m]33 


(14) 


where 
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Substituting  equations  (7),  (12)  and  (13)  into  (6),  the 
following  nonlinear  equation  of  motion  of  the  beam 
without  crack  is  obtained 

[m6  ] {? } + \Ke  (M)]  M = {pe } • (1 8) 

In  case  with  the  damping  force  j/J  j = ^Ce  » 
equation  (18)  becomes 


where 


Hc'] 

M- 
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>41[ 

+N( 

is  the  damping  matrix  of  element,  j^jis 
velocity  vector. 

2.2.  Beam  element  with  crack 

Considering  the  beam  element  with  crack,  stiffness  matrix 


of  the  element 

Kec({q 


K*r 


is  identified  as  [1 1],  [20] 
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where  \KC ] is  the  stiffness  matrix  of  weak  beam  element 

due  to  cracks,  and  the  above  matrices  can  be  formulated 
as: 
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with  components  in  equation  (21) 
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where  Jo,  Jc  , respectively,  are  the  moments  of  inertia  of 
the  beam  cross  section  for  Oz  axis  at  non-cracked 
positions  and  at  the  positions  of  cracks;  lc  = 1,3d  (d  is  the 
depth  of  the  crack),  le  is  the  length  of  the  element,  E is  the 
modulus  of  elasticity,  T|  is  the  distance  from  the  left  end 
of  beam  element  to  the  crack. 

Considering  that  the  lost  mass  due  to  cracks  is  little 
compared  to  the  overall  mass  of  the  element. 
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2.3.  Beam  element  on  elastic  foundation 

Stiffness  matrix  of  the  beam  element  on  elastic 


foundation  is  identified  by  [3]: 


K\ 


where 


Ke 


cracked  beam  element, 


(M)]-NM) 
■ M 


(22) 

correlates  with  the 
is  the  stiffness  matrix 


related  to  an  elastic  foundation. 

2.4.  Nodal  load  vector  element  beam  on  elastic 
foundation  under  moving  mass 

According  to  FEM  method,  when  a moving  load  is 
involved  in  the  working  of  the  system,  due  to  the  position 
change  property  of  the  load  over  time,  so  at  each  point  of 
time,  the  moving  load  acts  on  one  beam  element. 
Considering  the  beam  element  on  elastic  foundation 
subjected  to  the  moving  mass  m,  the  force  P(t)  acts  on  m 
(Fig.  3). 


Fig.  3.  Beam  element  on  elastic  foundation  under  moving 
mass. 

The  force  of  the  moving  mass  acting  on  the  beam  at  the 
coordinate  x = £,  = vt  is: 

l2',,  r , . 

(23) 

where  y(x,t)  is  element  deflection,  — is  the  absolute 

dr 

acceleration  in  the  y direction. 

After  taking  the  derivative  of  the  deflection  function,  the 
expression  (23)  is  rewritten  as 


/x  /x  d2y(x,t) 
R(x,t)  = P(t)-m — 

dt 1 


R^xA  - P(t)  — m 


( ~2  ^2  ^2  ^ 

dt2  - - 


dxdt  Qx2 


(24) 


The  acting  force  (24)  is  described  by  the  distributed  force 
p(x,t)  as: 

p(x,t)  = R(x,t)- S(x-vt).  (25) 

In  case,  beam  on  elastic  foundation: 
p(x,t ) = R(x,tf  5(x-  vt)-k$y(x),  (26) 

where  80  is  denotes  the  Dirac-Delta  function,  kO  - 

foundation  modulus. 

Therefore,  the  force  vector  is: 


(27) 


is  the  matrix  of  shape 


Fe)=\l‘[Nf  p(x,t)dx  = 

Jo  5^X  ~ £)[^]r  R (*’f  )*  - Jo  k0  [N  f y(x)dx, 

M [0]  [0] 

where,  [N]  = [0]  [0] 

[0]  [0]  \n^ 

functions  of  flexural  beam  element. 

Substituting  y (x)  = [TV]  j</ j into  Eq.  (24)  we  get 

(x,t)  = P(t)  -m([iV]{r}  + 2 v[N']{qe)  + [N"]{qe  j v2) 


R 


(28) 


(29) 


dN  d2N 

where  [N']=  — , [N"]=^4. 

fa  ax2 

Substituting  equations  (28)  into  (26)  and  paying  attention 
to  the  nature  of  the  Delta-Dirac  function,  equation  (27) 
may  be  rewritten  as 

where 

f (0}  = Jo  [Nf  s(x - (^)f  P (0.  (3°) 

\Mp\  = »jJ  S(x-Z)[N]T  [N]dx  = m[N(Z)]T  [lV(£)]  (31) 

| \cep\  = 2 vmjloeS(x  - ^)[N]T  [N']dx  = 2vm\_N  )], 

(32) 
le 

[ Kp ] = v2m\ S(x  - t)[N]T  [N"]dx  = mv2  [n(£)]T  [AT(f  )], 

0 

(33) 

[Kf]  = loko[N{x)f[N(x)]dx.  (34) 

Substituting  equation  (29)  into  equation  (19),  we  get  the 
equations  of  motion  governing  the  nonlinear  dynamic 
response  of  the  beam  element  on  elastic  foundation 
subjected  to  moving  mass 

([M']+[«;])(rj+([c<]+[c;])(4') 

2.5.  Governing  differential  equations  for  total  system 

Assembling  all  elements  matrices  and  nodal  force  vectors, 
the  governing  equations  of  motions  of  the  cracked  beam 
on  elastic  foundation  subjected  to  moving  mass  can  be 
derived  as 


(35) 


([M0]  + [Mp]){q}  + ([C0  ({<?})]  + [Cp]){q} 

+([%({«})] +M  M={p}’ 


(36) 
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where  [^o]  + [Mp]  = Z[M^]  + Z[Mp]’ 

e em 

[Q>  ({*})]  + [Cp ] = aR  [M0]  + Pr  [%  ({*})]  + Z [Cp\> 

em 

\j^bf  ] = l[Kbf\  ’ e *s  num^er  °f  normal  elements,  em 

e 

is  the  number  of  elements  directly  subjected  to  moving 
mass  and  aR,  Pr  are  the  Rayleigh  damping  coefficients. 
This  is  a nonlinear  differential  equation  system  with  time 
dependence  coefficient  that  can  be  solved  by  using 
Newmark’s  direct  integration  and  Newton-Raphson 
iteration  method.  A ANSYS  program  called 
CBF_Moving_Mass_2017  was  conducted  to  solve 
equation  (36).  The  code  of  the  calculation  program  is 
written  in  the  ANSYS  13.0  environment. 

III.  NUMERICAL  ANALYSIS 

Beam’s  length  L = 8m,  rectangular  cross  section  with 
width  b = 0.1m,  height  h = 0.2m;  one  end  is  in  pinned 
connection,  and  the  other  end  is  in  roller  connection. 
Beam  material  with  elastic  modulus  E = 2.1xlOn  N/m2, 
Poisson's  ratio  v = 0.3;  density  p = 7800  kg/m3  is  used. 
There  is  a V-shaped  open  crack  in  the  center  of  the  beam, 
and  the  crack’s  depth  d = 0.1  m.  Foundation  stiffness  ko  = 
2xl04  N/m3.  The  used  load  is  the  material  point  with  the 
mass  m = 1000kg,  moving  along  the  beam  with  the 
velocity  v = 36  km/h. 

With  the  established  program  established,  we  calculate 
for  02  cases:  Beam  with  crack  (basic  problem  - BP)  and 
without  crack  (comparison  problem  - CP)  to  see  more 
clearly  the  impact  on  the  dynamic  response  of  the  system 
when  the  cracks  appear.  The  response  results  of 
deflection  y,  acceleration  y , cutting  force  Qy  and  normal 

stress  gx  at  the  midpoint  of  the  beam  (point  A (4,0,0))  are 
shown  in  Table  2 and  Figures  4,  5,  6,  7.  Through  this 
results,  we  realize  that  with  cracked  beam  the  whole 
displacement,  acceleration  of  vertical  displacements  and 
normal  stresses  are  greater  than  the  beam  without  crack. 
This  problem  showed  the  dangers  of  crack  to  stiffness, 
stability  of  cracked  beam  on  elastic  foundation  under 
moving  loads. 

Table  1.  Summary  of  maximum  values  of  calculated 


quantities 


Quantities 

Ymax 

[cm] 

Ymax 

[m/s2] 

nmax 

[N] 

_max 

ax 

[N/m2] 

BP  (with 
crack) 

0.373 

20.994 

109.889 

3.316xl07 

CP  (without 
crack) 

0.191 

0.125 

119.421 

1.079xl07 

Differents 

1.95 

times 

167.95 

times 

0.92 

times 

3.07  times 

...•BP  CP 

0.002 
0.001 
I 0 

g 0 001 

r 0 002 
J 0003 
0 004 
0.005 

0,00  0.20  040  0.60  0.80 

Timet  | »| 

Fig.  4.  Response  of  deflection  y over  time  at  the  center 
cross  section  of  the  beam 
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Fig.  5.  Response  of  acceleration  y over  time  at  the  center 
cross  section  of  the  beam 


bp  CP 


100 

2 50 

f 0 
| 

If  *50 

3 *ioo 
*150 

0.00  020  0.40  0.60  080 

Time  t [sj 

Fig.  6.  Response  of  cutting  force  Qy  at  the  center  cross 
section  of  the  beam 
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Fig.  7.  Response  of  normal  stress  ox  at  the  center  cross 
section  of  the  beam 


The  results  show  that  the  effect  of  cracks  on  the  dynamic 
response  of  the  beam  is  significant.  For  cracked  beams, 
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vibration  of  the  beams  increases  after  the  mass  moves 
through  the  crack. 

3.1.  Effect  of  elastic  foundation  stiffness 

Studying  the  changes  in  maximum  values  of  the 
displacement,  internal  force  and  direct  stress  of  the  beam 
under  the  elastic  foundation  stiffness,  through  the 
stiffness  kO  of  the  spring  ranging  from  lxl04N/m3  to 
6xl04N/m3.  The  results  of  changes  in  maximum  values 
of  the  displacement,  internal  force  and  direct  stress  at  the 
center  cross  section  of  the  beam  are  shown  in  Table  3 and 
graphs  in  Figures  8 and  9. 

Table  3.  Summary  of  maximum  values  of  quantities 


calculated  based  on  ko 


ko  [N/m] 

lxlO4 

2xl04 

4xl04 

6xl04 

Ymax  [cm] 

0.459 

0.373 

0.279 

0.228 

Ymax  [m/s2] 

17.96 

3 

20.994 

21.762 

21.480 

Q“ax  [N] 

113.2 

20 

109.88 

9 

104.52 

5 

99.085 

<rfax  [N/m2] 
(xlO7) 

3.704 

3.316 

2.859 

2.590 

„ 0.55 
0.45 

c 0.35 

0 

1 0.25 
§ 0.15 

0 20000  40000  60000 

Foundation  stiffness  kO  [N/m] 


Fig.  8.  ymax  - ko  relation 
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Fig.  9.  G™  - ko  relation 

It  is  observed  that  when  the  foundation  stiffness 
increases,  the  maximum  values  of  displacement  and 
flexural  moment  decrease  due  to  the  increase  in  the 
system’s  overall  stiffness.  The  maximum  values  of 
displacement  and  flexural  moment  decrease  sharply  when 


ko  varies  from  lxl04N/m3  to  3xl04N/m3,  then  the 
decreasing  rate  shall  be  slower. 

3.2.  Effect  of  load  speed 

Surveying  the  problem  with  a load  speed  v changes  from 
lOm/s  (36km/h)  to  35m/s  (126km/h).  The  results  of  the 
variations  of  the  maximum  values  of  deflection, 
acceleration,  cutting  force  and  stress  at  the  midpoint  of 
the  beam  based  on  v are  shown  in  Table  4 and  graphs  in 
Figures  10  and  11. 

Table  4.  Summary  of  maximum  values  of  quantities 
calculated  based  on  v 


v [km/h]  36  90  126 


Ymax  [Cm] 

0.373 

0.290 

0.282 

Ymax  [m/s2] 

20.994 

23.102 

27.655 

Q“ax  [N] 

109.889 

90.804 

82.960 

a'™x  [N/m2] 

3.316 

2.890 

2.797 

(xlO7) 
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Fig.  10.  ymax  - v relation 
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Fig.  11.  a™ax  - v relation 

It  is  clear  that  when  the  moving  speed  of  the  load 
increases,  the  maximum  values  of  displacement,  internal 
force  and  stress  in  the  beam  decrease,  when  the  moving 
speed  of  the  load  varies  from  90  m/s  to  1 10  m/s,  the  direct 
stress  does  not  change  much,  then  decreases  sharply. 

3.3.  Effect  of  crack  location 

This  example  studies  the  changes  in  maximum  values  of 
displacement  and  internal  force  of  the  beam  according  to 
the  crack  location,  giving  the  cracks  located  far  from  the 


3.40E+07 


c „ 3.20E+07 

.1 " 

x J 3.00E+07 


CO 

CO 


aj  “ 2.80E+07 

■M 

1/1  2.60E+07 


www.ijaers.com 


Page  | 78 


International  Journal  of  Advanced  Engineering  Research  and  Science  (IJAERS) 
https: //dx.  doi.  org/1 0 . 221 61/ijaers.  4.9.14 


[Vol-4,  Issue-9,  Sep-  2017] 
ISSN:  2349-6495(9)  / 2456-1908(0) 


beam’s  ends  the  distances  L/4,  L/2,  3L/4.  Results  of 
extreme  values  of  the  responses  at  the  calculated  points 
are  shown  in  Tab.  5 and  graphs  in  Figures  12,  13,  14,  15. 
When  mass  moving  through  the  crack,  the  beam  vibrates 
amplitude  and  frequency,  which  shows  that  the 
stabilization  of  the  beam  oscillations  during  this  period 
decreased. 

Table  5.  Summary  of  maximum  values  of  quantities 


calculated  according  to  the  crack  location 


Crack  location 
(from  left  end) 

L/4 

L/2 

3L/4 

Ymax  [m] 

0.257 

0.373 

0.279 

ymax  l^l 

16.347 

20.994 

19.855 

Q“ax  [N] 

75.374 

109.889 

63.418 

o'™x  [N/m2] 

1.453 

3.316 

1.133 

(xlO7) 


Fig.  12.  Response  of  y according  to  the  crack  location 
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Fig.  13.  Response  of  y according  to  the  crack  location 
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Fig.  14.  Response  of  Qy  according  to  the  crack  location 


L/4 L/2  3 L/4 

x 5.00E+07 
| ^ 0.00E+00 

« c -5.00E+07 

X 

2 0.00  0.50 

i/> 

<u  — 

i:  Time  t [s] 

wo 

Fig.  15.  Response  ofax  according  to  the  crack  location 

Crack  location  changes  making  the  maximum  responses 
of  displacement,  stress  and  internal  force  in  the  beams 
change  significantly;  when  the  crack  is  in  the  center  of  the 
beam,  the  above  quantities  reach  the  maximum  values,  so, 
the  beam  is  most  dangerous  when  there  is  a crack 
appearing  in  this  position. 

IV.  CONCLUSION 

A conclusion  section  must  be  included  and  should 
indicate  clearly  the  advantages,  limitations,  and  possible 
applications  of  the  paper.  Although  a conclusion  may 
review  the  main  points  of  the  paper,  do  not  replicate  the 
abstract  as  the  conclusion.  A conclusion  might  elaborate 
on  the  importance  of  the  work  or  suggest  applications  and 
extensions. 

The  nonlinear  dynamics  analysis  of  cracked  beams 
resting  on  a Winkler  foundation  subjected  to  a moving 
mass  using  the  finite  element  method  has  been  presented. 
A two-node  beam  element  based  on  Euler-Bernoulli  beam 
theory,  taking  the  effect  of  crack  and  foundation  support, 
was  derived  and  employed  in  the  analysis.  The  dynamics 
response  of  the  beams,  including  the  time  histories  for 
deflection,  acceleration  and  normal  stress,  was  computed 
with  the  aid  of  Newmark  method.  The  effect  of  loading 
parameters,  foundation  stiffness  and  crack  location  on  the 
dynamic  response  of  the  beams  has  been  examined  and 
highlighted.  The  main  conclusions  can  be  summarized  as 
follows: 

The  beam  element  and  computer  code  developed  in  the 
present  work  are  accurate  in  evaluating  the  dynamic 
characteristics  of  cracked  beams  subjected  to  moving 
masses. 

The  elastic  foundation  plays  an  important  role  in  the 
dynamic  response  of  the  cracked  beams  under  a moving 
mass.  Both  the  dynamic  deflection  and  normal  stress  are 
significaly  decreased  by  the  increase  of  the  foundation 
stiffness. 

The  dynamic  response  of  the  cracked  beams,  as  in  case  of 
the  uncracked  beams,  is  governed  by  the  moving  mass 
speed.  With  the  moving  speeds  in  the  range  considered  in 
this  paper,  both  the  dynamic  deflection  and  normal  stress 
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decreased  when  increasing  the  moving  speed. 

The  maximum  dynamic  deflection  and  normal  stress  are 
significantly  influenced  by  the  crack  location.  The 
deflection  and  normal  tress  attain  the  largest  values  when 
the  crack  is  located  at  the  midpoint  of  the  beam.  Thus, 
from  an  engineering  point  of  view,  the  midpoint  crack  is 
the  most  dangerous  one. 

The  results  obtained  in  this  paper  help  to  select 
appropriate  parameters,  the  solution  for  structural 
reinforcement  cracked  beam  on  elastic  foundation  under 
moving  load  and  applications  in  transportation  techniques 
such  as  the  train  rails. 
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